NPS  ARCHIVE 
1962 
STECKBECK,  F. 


_nB§l 


•  ••■  ■■■''■■ 

191 

Mil   :  .  '  '  . 


A   STATISTICAL  METHOD   FOR   FORECASTING  THE 
MOVEMENT  OF  NORTH  AMERICAN  ANTICYCLONES 

FRANCIS  J.  STECKBECK 

and 

ALBERT  H.  MANHARD,  JR. 


'V- 


W** 


JPBHS 

■  ■■■  ■'"■'■- '-■■'■ 
■'■■''■'■■■. 

I '  ■■  H 

',:".■.'•'. 

am 

-.■■■'"•■. 


mm 

tflSffit 

'•■:'•■ 

'■■•": 

;(•'.•■.■•■■'. 


LIBRARY 

U.S.  NAVAL  POSTGRADUATE  SCHOOL 

MONTERFY.  CALIFORNIA 


A  STATISTICAL  METHOD  FOR  FORECASTING 
THE  MOVEMENT  OF 
NORTH  AMERICAN  ANTICYCLONES 


•sV  -/V  Vf  -,'c  •&■ 

Francis  J.  Steckbeck 

and 
Albert  H.  Manhard  Jr 


A  STATISTICAL  METHOD  FOR  FORECASTING 
THE  MOVEMENT  OF 
NORTH  AMERICAN  ANTICYCLONES 


by 
Francis  J.  Steckbeck 

9 

Lieutenant  Commander,  U.S.  Navy 
and 
Albert  H.  Manhard  Jr. 
Captain,  U.S.  Marine  Corps 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

METEOROLOGY 


United  States  Naval  Postgraduate  School 
Monterey,  California 

19  6  2 


^ 

-> 


^ecWJjp 


-e^fc- 


LIBRARY 

U.S.  NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CALIFORNIA 


A  STATISTICAL  METHOD  FOR  FORECASTING 
THE  MOVEMENT  OF 
NORTH  AMERICAN  ANTICYCLONES 
by 
Francis  J.  Steckbeck 

and 
Albert  H.  Manhard  Jr, 

This  work  is  accepted  as  fulfilling 

the  thesis  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

METEOROLOGY 

from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 

A  statistical  prediction  method  to  forecast  the  movement  of  North 
American  anticyclones  is  developed  utilizing  multiple  linear  regression 
analysis  performed  by  the  CDC  1604  electronic  computer.   The  three  pre- 
dictands  sought  are  the  24-hour  change  in  central  pressure  and  the  east- 
ward and  northward  displacements.   Three  multiple  linear  regression 
equations  are  derived. 

The  input  data  consist  of  observations  of  150  anticyclones.   A 
moving  coordinate  system  is  employed:   the  observed  and  derived  pre- 
dictor information  is  measured  relative  to  the  system  center  rather  than 
a  fixed  geographical  position. 

Readily  measured  meteorological  parameters  are  selected  for  the 
input  data.   These  include  point  values  of  surface  pressure,  surface 
temperature,  500-mb  height  and  average  temperature  of  the  layer  from  the 
surface  to  the  500-mb  level.   The  24-hour  change  of  these  variables  is 
also  included. 

The  regression  equations  are  applied  to  50  independent  test  cases. 
These  statistical  forecasts  incur  a  vector  error  of  254  nautical  miles 
for  position,  and  a  mean  absolute  error  of  3.5  mb  for  pressure  change 
at  the  system  center. 

The  report  is  modeled  after  a  similar  study  of  cyclones  conducted 
by  Fredrick  P.  Ostby  and  Keith  W.  Veigas  [i960]. 

The  authors  are  indebted  to  Professor  Frank  L.  Martin  for  his  advice 
and  guidance  during  the  progress  of  this  study.   Acknowledgment  is  also 
due  Professor  J.  R.  Borsting  for  his  invaluable  assistance  in  the 
statistical  analysis  of  the  investigation. 
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1.    Introduction, 

Recent  advances  in  statistical  theory  combined  with  the  per- 
fection of  the  modern  electronic  computer  make  statistical  methods 
useful  in  nearly  every  area  of  scientific  endeavor.   Such  machine 
analyses  are  readily  adaptable  to  the  collection,  investigation  and 
interpretation  of  quantitative  data  in  a  myriad  of  fields.   Logical 
statistical  analysis  is  especially  useful  to  the  meteorologist. 

The  meteorologist  collects  data  which  are  relatively  complex  and 
which  are  the  result  of  numerous  factors.   Thus,  occurrences  of 
meteorological  phenomena  are  not  independent.   Statistical  predictions 
are  only  average  values  based  upon  the  particular  period  examined.   We 
never  do  know  the  true  probability  of  a  meteorological  phenomenon;  we 
can  only  form  an  estimate  of  it  based  on  past  experience.   Essentially 
accurate  conclusions,  of  course,  may  be  drawn  from  rough  estimates,  but 
in  their  use  the  element  of  danger  is  correspondingly  greater.   Sound 
judgment  must  be  constantly  invoked  to  guard  against  false  conclusions 
being  drawn. 

The  physicist  and  the  biologist  utilize  laboratory  techniques  to 
deal  with  problems  of  multiple  relationships.   Under  laboratory  condi- 
tions all  the  variables  except  the  one  under  consideration  can  be  held 
constant  while  the  latter 's  effect  upon  the  dependent  variable  can  then 
be  determined.   For  the  problems  of  the  meteorologist,  however,  such 
laboratory  controls  cannot  be  imposed.   Rainfall,  cloud  cover,  tempera- 
ture and  pressure  vary  constantly  and  only  their  combined  effect  can 
be  noted.   Hence,  the  meteorologist  must  frequently  apply  methods  of 
multivariate  statistical  analysis  in  order  to  handle  serially- 
correlated  data,, 
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Statistics  may  be  collected  on  most  observations,  but  the  use  of 
statistical  methods  in  these  studies  does  not  always  lead  to  results  of 
universal  validity.   At  best  the  statistical  may  be  only  one  of  many 
approaches  to  the  explanation  of  physical  phenomena.   Furthermore, 
statistical  inferences  must  always  be  scrutinized  with  care.   On  the 
other  hand,  the  successful  application  of  statistics  and  statistical 
methods  to  meteorological  problems  may  be  sufficiently  emphasized,  at 
this  point,  by  calling  attention  to  only  a  few  of  the  significant 
statistical  studies  in  the  field.   Useful  contributions  have  been  made 
by  Charney  [1954],  Miller  [1958],  Ostby  and  Veigas  [1960],  and  Hull  [1961]. 

The  aim  of  this  study  is  to  utilize  statistical  methods  to  develop 
multiple  linear  regression  equations  for  predicting  the  behavior  of 
migratory  anticyclones  in  North  America.   Such  equations  would  be  use- 
ful tools  to  both  the  synoptic  and  numerical  weather  forecaster. 
2.    Background. 

The  development  of  anticyclones,  from  the  standpoint  of  the  mass 
and  temperature  fields,  is  the  reverse  of  that  of  cyclonic  systems. 
High-level  anticyclogenesis  appears  to  be  associated  with  the  accumula- 
tion of  mass  in  the  lower  stratosphere,  accompanied  by  cooling.   This 
cooling  can  be  advective  or  due  to  ascent  resulting  from  horizontal 
convergence  in  the  upper  troposphere.   For  a  complete  discussion  of  the 
theories  dealing  with  the.  origin  and  structure  of  anticyclones,  the 
reader  is  referred  to  an  article  by  Wexler  [1951]. 

Harrison  [1947],  investigating  the  steering  of  anticyclones,  found 
that  about  687,  of  the  North  American  anticyclones  were  steered  by  the 
flow  at  700  mb.   Longley's  [19'v7]  results  showed  that  about  757«  were 
steered  by  this  same  current.   However,  if  either  of  the  above 
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techniques  are  to  be  applied  to  predict  the  movement  of  a  surface  anti- 
cyclone, first  the  700-mb  flow  pattern  must  be  correctly  forecast.   This 
of  course  reduces  the  value  of  the  methods  to  the  forecaster. 

George  and  associates  [1960]  have  devised  a  semi-objective  technique 
for  forecasting  the  movement  of  anticyclones  for  24-hour  and  48-hour 
periods.   This  solution  uses  the  isotherm  ribbon  at  850  mb  for  determin- 
ing the  direction  of  movement,  and  the  first  definite  wind  maximum,  "at 
700  mb  transverse  to  the  current  flow  adjacent  to  the  surface  high 
center,"  to  determine  the  speed  of  movement.   An  independent  data  check 
by  George  et  al.  showed  the  average  error  in  position  to  be  265  miles  for 
the  24-hour  forecast.   However,  when  utilizing  this  method,  the  flow 
pattern  and  its  relation  to  the  surface  anticyclone  remain  somewhat  sub- 
jective. 

Successful  objective  forecast  schemes  are  of  considerable  value  in 
meteorology.   In  designing  such  a  forecast  procedure,  one  should,  if 
possible,  include  as  quantitative  predictors  those  parameters  considered 
important  on  a  theoretical  or  semi-theoretical  basis.   However,  loose 
reasoning  may  result  if  one  attempts  to  assign  cause  and  effect  based 
upon  statistical  investigations.   In  this  study  no  statement  concerning 
cause  and  effect  will  be  made.   Rather,  the  results  will  be  judged  on 
the  performance  of  the  derived  regression  equations  themselves. 
3.    Data  sources. 

The  anticyclonic  cases  used  in  this  study  were  taken  from  the  Daily 
Series  Synoptic  Weather  Maps  analyzed  by  the  National  Weather  Analysis 
Center  for  both  the  surface  and  500-mb  levels.   These  charts  use  a  Polar 
Stereographic  Projection  true  at  latitude  60N„   They  were  utilized  for 
compiling  both  the  dependent  and  independent  test  data.   The  sea-level 
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maps  were  analyzed  for  observations  taken  at  1230Z,  while  the  500-mb 
observations  were  taken  at  1500Z. 

Anticyclones  for  the  dependent  sample  were  selected  from  a  geo- 
graphic area  extending  from  longitude  60  to  HOW,  between  latitudes  30 
and  60N.   This  area  is  outlined  in  figure  1.   The  period  covered 
extended  from  October  through  March  for  the  years  1953  through  1957 
(1955  excluded).   The  final  regression  equations  developed  in  this  study 
should  be  applied  only  from  October  through  March  for  those  anticyclones 
which  lie  in  the  area  defined  above,  henceforth  called  the  forecast  area. 

In  order  for  an  anticyclone  to  be  included  in  the  sample,  it  had  to 
be  characterized  by  at  least  one  closed  isobar  (5-mb  interval),   A 
further  requirement  was  that  the  closed  system  had  to  remain  within  the 
forecast  area  during  the  succeeding  24  hours.   If  the  system  met  both  of 
the  above  requirements  for  an  additional  24-hour  period,  it  was  included 
again  in  the  sample  as  another  case.   In  order  to  eliminate  any  possible 
bias,  no  single  system  was  included  a  third  time.   These  precautions  are 
similar  to  those  taken  by  Ostby  and  Veigas  [1960]  in  their  study  of 
cyclones. 

Based  upon  the  above  criteria,  150  cases  were  selected  for  the 
dependent  sample  data.   The  frequency  distribution  for  this  data  by  month 
and  year  is  as  shown  in  table  1.   The  number  of  cases  examined  varied 
from  as  many  as  11  in  February  of  1953  to  as  few  as  one  in  February  of 
1957.   The  number  of  cases  was  a  maximum  during  March  (an  average  of  8.2) 
and  a  minimum  in  February  (an  average  of  6.0). 

The  year  1955  was  not  included  in  the  dependent  sample  in  order  that 
the  data  for  this  period  might  be  retained  for  use  in  the  independent 
test. 
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Oct 

Nov 

Dec 

Jan 

Feb 

Mar 

Total 

1953 

_ 

9 

11 

9 

29 

1953-54 

-- 

9 

7 

8 

5 

10 

39 

1954-55 

5 

10 

10 

-- 

-- 

-- 

25 

1955-56 

-- 

-- 

-- 

6 

7 

9 

22 

1956-57 

8 

6 

8 

7 

1 

5 

35 

Total 

13 

25 

25 

30 

24 

33 

150 

Table  1.   Dependent  sample  frequency  distribution  by  months  and  years. 


4.    Selection  of  possible  predictors. 

The  development  of  the  modern  electronic  computer  has  made  it 
possible  to  analyze  statistically  a  large  number  of  possible  predictors 
with  great  speed.   A  screening  program  adapted  to  the  CDC  1604  electronic 
computer  is  one  method  of  carrying  out  the  preliminary  analysis  of  the 
reduction  of  variance  necessary  for  the  selection  of  predictors.   This 
program  is  a  modification  of  one  originally  written  by  M.  A,  Efroymson 
[1960]  of  the  Esso  Research  and  Engineering  Company „ 

When  utilizing  such  a  technique,  the  selection  of  predictors  must 
be  purely  objective.   Consequently,  it  is  imperative  that  the  user 
specify  precisely  the  initial  set  of  possible  predictors.   The  considera- 
tion of  known  synoptic  or  dynamic  influences  should  be  a  factor  in  the 
selection  process.   Dynamic  principles  must  be  considered  if  any  signifi- 
cant degree  of  success  is  to  be  expected. 

A  grid  system  capable  of  specifying  point  values  of  many  atmospheric 
variables  over  the  forecast  area  was  devised,   A  moving  coordinate  system 
was  selected  similar  to  that  of  Ostby  and  Veigas  [i960].   The  observed 
and  derived  predictor  information  is  measured  relative  to  the  anticyclone 
center  rather  than  a  fixed  geographical  position.   The  basic  grid  is 
shdwn  in  figure  2.   It  should  be  noted  that  this  grid  differs  from  that 
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Figure  2.   The  basic  grid*   The  y-axis  is  aligned  along  a  meridian  through 
the  sea-level  anticyclone. 
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of  Ostby  and  Veigas  in  being  square  rather  than  determined  by  a  latitude- 
longitude  net.   This  square  grid  was  constructed  by  first  locating  grid 
point  54  along  latitude  45N.   All  remaining  grid  points  are  then  located 
at  5-degree  latitude  intervals  both  in  the  x  and  y  directions.   This 
5-degree  interval  is  actually  300  nautical  miles  measured  at  45N.   The 
basic  grid  was  utilized  for  determining  all  parameters.   Grid  point  54 
is  located  over  the  system  center  with  the  y-axis  aligned  along  the 
meridian.   Smaller  sets  of  grid  points  have  been  used  for  recording 
certain  parameters,  but  these  represent  merely  sub-sets  of  the  basic  grid 
and  utilize  the  same  numbering  system  for  the  various  grid  points. 

Features  of  the  present  and  past  24-hour  pressure  and  temperature 
fields  are  embodied  to  furnish  the  majority  of  possible  predictors.   The 
sea-level  chart  and  the  500-mb  contour  chart  were  employed  to  obtain  the 
various  possible  pressure-height  predictors.   The  procedures  employed  are 
listed  below  under  the  headings  (a)  sea-level  predictors  and  (b)  upper -air 
predictors. 

(a)   Sea-level  predictors.   Point  values  of  the  sea-level  pressure  field 
were  recorded  and  included  as  possible  predictors.   The  grid  overlay  was 
positioned  on  the  surface  chart  with  grid  point  54  located  over  the 
system  center.   Values  of  sea-level  pressure,  read  to  the  nearest  whole 
millibar,  were  then  recorded  at  the  54  grid  points  as  indicated  in 
figure  3.   These  values  then  specified  the  sea-level  pressure  field  around 
the  system  center  at  observation  time  (t^) . 

With  the  grid  in  the  same  relative  position,  point  values  of  surface 
pressure  24  hours  prior  to  the  observation  time  (t_-24)  were  recorded  at 
the  13  shaded  grid  points,  as  also  indicated  in  figure  3.   These  values 
were  then  used  to  specify  13  values  of  the  24-hour  pressure-change  field, 

-8- 


> 

r 

o 

12 

o 

22 

• 

52 

o 

42 

I     o 

52      62 

• 

72 

O 

82 

O 

92 

o 

15 

o 

25 

o 

55 

• 

45 

O        • 

55      65 

O 

75 

O 

85 

O 

91 

o 

A 

-  O 

— o — 

44 

-• — o— 

54      64 

— C— 

74 

— #— 

64 

—O- 

94 

14 

24 

54 

o 

15 

o 

25 

o 

55 

• 

45 

C 
5! 

5     65 

O 

75 

O 

85 

O 

95 

o 

16 

o 

26 

# 

56 

o 

46 

4 

5< 

►     O 

5      66 

76 

O 

86 

c 

96 

o 

17 

o 

27 

O 
57 

o 

47 

C 

5 

3      O 

7      67 

O 

77 

O 

r>,~7 

O 

97 

• 

n 

at  tQ-24 

Figure  5.   3rid  coints  used  tc  obtain  sep-level  pressures. 
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which  were  also  included  in  the  list  of  possible  predictors. 

The  latitude,  longitude  and  intensity  of  the  anticyclone  at  t_  were 
also  recorded  for  use  as  possible  predictors.   The  intensity  was  measured 
by  averaging  the  surface  pressure  at  grid  points  43,  45,  63  and  65,  then 
subtracting  this  average  from  the  surface  pressure  at  grid  point  54, 
(b)   Upper  air  predictors.   The  values  of  the  500-mb  predictors  at  t 
were  obtained  from  the  500-mb  chart  by  positioning  grid  point  54  of  the 
appropriate  grid  over  the  surface  anticyclone  center  at  t  „   Values  of 
500-mb  heights  were  then  recorded  at  27  grid  points  as  indicated  in 
figure  4.   These  27  values  were  also  included  in  the  list  of  possible 
predictors. 

With  the  grid  in  the  same  relative  position,  point  values  of  500-mb 
height  24  hours  prior  to  the  observation  time  were  recorded  at  the  same 
27  grid  points.   Recorded  values  were  then  used  to  determine  the  24-hour 
500-mb  height  change  field.   These  values  were  then  included  as  possible 
predictors. 

As  previously  stated,  high-level  anticyclogenesis  appears  to  be 
associated  with  the  accumulation  of  mass  in  the  upper  troposphere  and 
lower  stratosphere,  together  with  associated  cooling.   Using  the  thermal 
pattern  implied  by  a  grid  plot  of  temperature  averaged  in  the  vertical, 
one  may  use  the  dynamical  concepts  of  Sutcliffe  [1947]  to  deduce  at  least 
qualitative  approximations  to  the  upper-level  convergence,  and  hence  the 
evolution  of  the  synoptic  flow  pattern  of  middle  and  upper  layers  of  the 
troposphere.   According  to  Sutcliffe,  baroclinic  development  and  movement 
of  a  sea-level  anticyclone  may  be  dynamically  related  to  its  position 
relative  to  the  1000-500  mb  mean  temperature  field*   In  order  to  arrive 
at  a  representative  mean-temperature  field,  the  surface  and  500-mb 
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temperatures  were  recorded  at  19  grid  points,  as  indicated  in  figure  5„ 
From  these  values  the  mean  temperatures  of  the  layer  were  calculated  and 
also  included  in  the  list  of  possible  predictors. 

It  should  be  noted  that  we  are  seeking  linear  regressions,  whereas 
Sutcliffe's  computations  of  advection  of  vorticity  or  of  thermal  vorticity 
involve  non-linear  fields.   Any  linear  regressions  obtained  here  which 
are  consistent  with  his  philosophy  cannot  be  expected  to  give  the  degree 
of  significance  that  would  exist  using  appropriate  non-linear  regressions. 

A  total  of  145  possible  predictors  plus  three  predictands  were 
recorded.   The  possible  predictors  are  summarized  in  table  2. 


Number  of 
Predictor  Symbol  predictors 

Sea-level  pressure  p  54 

24-hour  surface  pressure  change  AP  13 

Surface  temperature  T  1 

24r-hour  surface  temperature  change  AT  1 

Location  of  surface  anticyclone  L,  ?<  2 

Anticyclone  intensity  I  1 

500-mb  height  z  27 

24-hour  500-mb  height  change  Az  27 

Mean  temperature,  surface  to  500  mb  T  19 

145 


Table  2.   Possible  predictors  included  in  the  dependent  data. 

The  data  were  checked  and  verified,  then  punched  onto  IBM  cards. 
These  cards  were  also  checked  and  corrections  inserted  where  discrep- 
ancies existed.   The  data  were  then  transferred  to  magnetic  tape,  where 
they  were  stored  in  a  form  readily  accessible  for  repeated  machine 
analysis. 
5.   Reduction  of  data. 

T.e  reduction  of  the  145  possible  predictors  was  accomplished  by 
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employing  a  stepwise  regression  procedure.   The  regression  equations  were 
developed  by  adding  one  independent  variable  at  a  time  at  each  step  of 
the  calculation.   In  each  case,  the  variable  added  was  that  one  which 
made  the  r.ost  significant  improvement  in  "goodness  of  fit,"  based  upon  an 
F  test.   An  important  feature  of  the  regression  procedure  is  that  a 
variable  may  appear  to  be  significant  in  the  early  stages  (and  thus  enter 
the  linear  statement  of  the  problem),  but  after  several  other  variables 
are  added  to  the  regression  equation  prove  to  be  insignificant.   The 
insignificant  variable  is  then  removed  from  the  equation  before  another 
variable  is  added.   This  may  happen,  for  example,  when  the  variable 
initially  added  is  approximately  a  linear  combination  of  several  other 
independent  variables  which  have  subsequently  been  introduced  into  the 
regression  equation.   Thus,  only  significant  variables  are  retained  for 
inclusion  in  the  final  equation.   Additional  variables  are  then  added  in 
the  same  manner.   This  procedure  is  repeated  until  each  remaining 
variable  examined  fails  to  explain  a  specified  additional  percentage  of 
the  remaining  variance  of  the  predictand. 

The  F  level  for  both  entering  and  removing  variables  was  specified 
at  3.92.   For  a  description  of  the  significance  of  this  value  of  the 
F  level,  see  section  7.   The  program  then  computed  the  predicted  values 
and  the  deviations  from  the  actual  values,  as  well  as  the  RMS  error. 

When  we  are  dealing  with  multivariate  analysis,  estimation  of  com- 
ponent variances  is  sometimes  useful.   However,  in  the  case  of  several 
variables  added,  this  was  not  possible  due  to  the  non- independence  of  the 
predictors.   This  aspect  is  also  dealt  with  in  section  7. 
6.   The  form  of  the  forecast  equations. 

For  forecasts  of  this  kind  to  be  useful,  it  is  necessary  that  a 
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change  in  the  variable  to  be  predicted  be  attended  by  some  corresponding 
change  in  the  other  variables.   Based  upon  section  5,  we  can  express  the 
value  of  the  dependent  variable,  Y,  in  terms  of  the  several  independent 
variables,  X..   The  multiple  regression  equation  for  obtaining  the  pre- 
dicted value  of  Y  may  be  written 

Y»   =  aQ  +  a1X]_   +  a£X2   +   .  .  .   +  aj^  (1) 

where  Y'  is  the  predicted  value  of  Y  for  each  set  of  values  for  X,  the 

a.  are  constants  derived  from  the  input  data  and  the  X . ' s  are  the  pre- 
i  i 

dictors  selected  by  the  screening  process. 

In  the  final  derived  regression  equations,  the  symbol  Y8  denotes  in 
turn: 

s  Predicted  24-hour  south-north  displacement  of  the  system 
center  in  nautical  miles.   Positive  (negative)  values  indi- 
cate northward  (southward)  displacement, 
E'  -     Predicted  24-hour  west-east  displacement  of  the  system  center 
in  nautical  miles.   Positive  (negative)  values  indicate  east- 
ward (westward)  displacement. 
=  Predicted  24-hour  change  in  central  pressure  in  millibars. 
Positive  (negative)  values  indicate  pressure  increases 
(decreases) , 
The  final  derived  regression  equations  are  summarized  in  table  7,  along 
with  some  intermediate  truncated  forms  of  these  equations. 

In  the  prediction  equations  each  predictor  is  identified  by  a 
symbol  followed  by  a  subscript.   In  each  case  the  subscript  refers  to  the 
location  of  a  grid  point  as  specified  in  figure  2*   The  units  appropriate 
to  each  predictor  are  described  in  Appendix  B, 

The  percentage  of  total  variance  of  the  predictand  explained  by  the 
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selection  of  each  predictor  is  summarized  in  table  3.   The  predictors  are 
listed  in  the  order  of  their  selection,  which  in  turn  is  a  result  of  the 
computer  output.   Generally,  the  first  two  or  three  predictors  explain  a 
substantial  part  of  the  variance  of  the  dependent  variable.   Each 
additional  predictor  then  accounts  for  a  lesser  amount. 
7.    Statistical  analysis  and  results. 

The  BIMD  09  of  the  BIMD  series  of  statistical  programs  was  employed 
by  the  CDC  1604  electronic  computer  to  perform  this  regression  analysis. 
The  BIMD  09  program  of  the  U.  S.  Naval  Postgraduate  School  Computer  Center 
is  a  modification  of  the  BIMD  09  program  of  the  University  of  California 
at  Los  Angeles.   The  latter  program  is  in  turn  a  modification  of  an 
original  program  developed  by  M.  A,  Efroymson  [1960]  of  the  Esso  Research 
and  Engineering  Company.   A  brief  description  of  the  program  has  been 
given  in  section  5  of  this  report. 

Miller  [1958]  has  developed  a  screening  procedure  in  which  the 
determination  of  the  F  level  for  significance  of  the  contribution  of  a 
new  independent  variable  depends  upon  the  number  of  possible  predictors 
and  the  number  of  selected  predictors.   This  particular  program  is  in 
machine  language  for  the  IBM  704  computer  and  is  presently  not  available 
for  the  CDC  1604.   Statistical  analysis  of  meteorological  parameters 
utilizing  Miller's  procedure  has  shown  promising  results  in  connection 
with  movement  and  deepening  of  cyclones. 

An  important  phase  of  the  BIMD  09  program  occurs  in  the  selection  of 
critical  F  levels  to  enter  and  remove  variables.   Once  selected,  these 
levels  remain  constant  and  do  not  change  each  time  that  a  variable  is 
selected  (as  it  does  in  Miller's  procedure).   Therefore,  although  the 
criterion  for  selection  is  initially  less  stringent  than  Miller's,  the 
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statistical  significance  of  each  variable  selected  is,  in  general, 
clearly  evident.   The  F  level  input  necessary  for  a  critical  F  limit  is 
a  function  of  the  number  of  variables  to  be  selected,  which  is  an  unknown 
until  the  regression  analysis  is  completed. 

Initially  the  F  levels  for  entering  and  removing  variables  was 
specified  at  4.5  and  0.22  respectively.   The  resulting  three  regression 
equations  contained  six,  nine  and  nine  variables  respectively  for  the 
south-north,  west-east  and  pressure-change  equations. 

With  the  above  information  then  available,  a  table  of  F  values  was 
entered  for  a  critical  F  corresponding  to  a  confidence  level  of  0.95  and 
with  one  and  140  degrees  of  freedom.   This  resulted  in  a  critical  F  level 
of  3.92  which  was  then  used  in  the  program  for  both  entering  and  removing 
variables. 

With  the  final  computer  print-out  of  the  regression  equations  and 
the  statistical  data  completed,  the  data  were  re-examined  and  the  new 
cumulative  percent  reductions  of  variances  were  determined,  and  found  to 
be  consistent  with  the  initial  reductions.   These  statistics  and  others, 
as  described  below,  are  listed  in  table  3. 

The  overall  F  levels  for  the  final  equations  were  desired.   These 
values  would  indicate  the  significance  of  the  predictive  power  of  the 
final  equations. 

Fisher  [1948]  describes  the  F  test  for  significance  of  the  multiple 
correlation  coefficient  as  the  ratio  of  the  mean  squares  explained  to  the 
mean  squares  residual  (as  also  used  by  Miller  [1958]).   Anderson  [i960] 
and  Panofsky  and  Brier  [1958]  give  a  similar  test  for  the  sample  multiple 
correlation  coefficient  with  independent  variables  from  a  multivariate 
normal  distribution,  as  an  F  test  consisting  of  a  ratio  of  two  %  - 
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distributed  variables  with  k  and  n-k-1  degrees  of  freedom.   Here  n  is  the 
sample  size  and  k  is  the  number  of  variables  in  the  regression  equation. 
The  form  of  the  F  test  used  to  compute  significance  of  kth  regression 
equation  was  a  modification  of  a  formula  of  Anderson's  [1960,  p.  89],  and 
is  given  in  equation  (2)  below: 


F, 


7? 


2. 


"?, 


t-A-i 


(k,  n-k-V  1-7?*  ft 

Here  R  is  the  multiple  linear  correlation  coefficient.   With  a  multiple 
linear  correlation  R,  the  reduced  standard  error,  S  ,  relative  to  the 
total  standard  error,  s  ,  is 

y 


(2) 


s     - 


-s  \ 


1  - 
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or,  solving  for  R  , 
Therefore,  we  have 
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.    y 


=  percent  reduction  of  variance.    (4) 
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which  was  the  form  of  the  equation  used  to  compute  the  significance  test 
for  the  final  regression  equations.   For  each  variable  added  by  the 
screening  procedure,  the  computer  print-out  included  the  new  value  of  S  . 

The  cumulative  percent  reduction  of  variance  was  then  computed  for 
each  step  of  the  regression  using  equation  (4).   Thus,  by  subtraction  of 
the  previous   cumulative  percent  reduction  of  variance,  one  obtains  the 
percent  reduction  of  variance  attributable  to  the  variable  just  entered. 
Note  that  in  table  3  no  P.R.  is  attributed  to  the  first  two  variables 
individually,  but  only  in  combination.   This  was  due  to  the  fact  that  the 
variable  z,_  was  entered  at  the  first  step  of  the  regression,  but  was 
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removed  after  the  third  step  was  completed.   The  F  level,  for  that 
variable  individually,  had  dropped  below  the  0.95  significance  level. 
This  made  it  extremely  difficult  to  compute  the  percent  reduction  of 
variance  for  the  first  variable  in  the  final  regression  equation.   A 
different  program  would  have  been  necessary  to  compute  this  percent 
reduction  of  variance,  but  this  was  not  done  here. 

The  multiple  linear  correlation  coefficient  was  also  computed  for 
the  final  equations.   It  is  merely  the  square  root  of  the  final  cumula- 
tive percent  reduction  of  variance. 

A  test  was  undertaken  to  determine  whether  the  coefficients  in  the 
final  regression  equations,  taken  individually,  were  significantly 
different  from  zero.   This  could  be  accomplished  easily  since  the  com- 
puter print-out  included  the  standard  error  of  the  regression  co- 
efficients.  Dividing  the  coefficient  by  its  standard  error  yields  the 
"t"  statistic  which  was  used  to  determine  the  "t"  significance  level 
listed  in  table  3.   Note  that  all  coefficients  but  the  last  have  signifi- 
cance levels  greater  than  0.99.   It  should  be  noted  that  the  values  to 
the  left  of  the  double  vertical  line  in  table  3  were  values  computed  at 
each  stage  of  the  regression  while  the  "t"  values  to  the  right  of  this 
line  apply  only  to  the  final  regression  equation. 

These  predictive  equations  were  then  tested  on  an  independent  sample 
of  50  anticyclones.   The  predicted  values,  actual  values  and  errors  are 
listed  in  table  6.   Various  simple  statistical  parameters  were  calculated 
by  the  computer  as  the  50  forecasts  were  generated.   From  these  simple 
statistical  parameters,  equation  (4)  was  used  to  compute  the  coefficient 
of  determination  (P,R.  of  variance)  and  the  multiple  linear  correlation 
coefficient.   These  values  are  listed  in  table  4. 
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S-N  W-E  Change  in 

Displacement     Displacement     central  pressure 


s 
y 

221.12  nm 

195.38  nm 

4.19  mb 

s 

y 

271.64  nm 

217.70  nm 

5.34  mb 

R 

.5809 

.4411 

.6193 

R2 

.3374 

.1946 

.3835 

P.R.  variance 

33 . 74% 

19.467. 

38.35% 

Mean  of  actual 

-111  nm 

383  nm 

0.5  mb 

Range  of  actual 

-1000  to  400 

-55  to  750 

-10  to  15 

Range  of  predicted 

-462  to  349 

-112  to  760 

-7.5  to  9. 

0 

Range  of  error 

-851  to  442 

-375  to  451 

-6.4  to  9, 

,1 

Mean  of  errors 

0.5  nm 

14.5  nm 

0.4  mb 

Mean  absolute  error 

164.5  nm 

158.4  nm 

3.5  mb 

Table  4.   Summary  of  statistical  parameters  of  independent  data. 

A  "t"  test  was  made  to  determine  whether  both  samples  were  from  the 
same  infinite  population  using  a  significance  level  of  0.95.   This  test 
is  described  in  Spiegel  [1961],   As  stated  by  Hoel  [1956],  the  validity 
of  this  test  is  dependent  upon  whether  the  two  samples  have  the  same 
standard  errors.   Both  propositions  verified  at  the  significance  level 
of  0.95. 

The  errors  in  movement  were  plotted  graphically  in  order  to  obtain 
a  mean  vector  error  magnitude  and  thus  a  visual  representation  of  the 
deviations  from  the  actual  movement.   The  magnitude  of  the  mean  vector 
error  was  254  nautical  miles.   This  error  distribution  is  shown  in 
figure  6. 

Note  that  for  the  dependent  sample  the  percent  reduction  of  variance 
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555/852  nm 


'tic^l  miles 


Firure  6, 


^rror  distribution  of  50  independent  predictions  of  24-hour 
anticyclone  displece  ent.   Cripin  represents  actual  locations, 
The  predicted  locations  rre  then  represented  by  small  circles 
plotted  relr  i- c  to  the  -ctual  locations.   The  5C>  °?  "nc*  95 
percentiles  are  shown, 
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for  all  three  equations  is  generally  507o.   This  reduction  compares  favor- 
ably with  Ostby  and  Veigas'  results  [1960]  for  the  more  predictable 
cyclonic  systems  in  a  considerably  smaller  forecast  area. 

The  independent  test  shows  a  small  loss  of  percent  reduction  of 
variance  except  for  the  west-east  equation.   This  equation  explained  307o 
less  of  the  percent  reduction  of  variance  when  tested  on  the  independent 
sample.   This  could  be  accounted  for  by  the  possibility  that  the 
variables  are  not  normally  distributed.   In  addition,  the  relationship 
between  the  predictand  and  the  predictors  may  not  be  linear.   A  better 
fit  could  perhaps  be  expected  from  a  non-linear  combination  of  these  or 
other  parameters.   Such  transformations  can  be  programmed  into  the 
BIMD  09  program,  although  this  was  not  within  the  scope  of  the  present 
work. 
8.   A  discussion  of  the  errors. 

Two  types  of  errors  are  considered  in  this  section:   those  pertaining 
to  data  acquisition;  and  prediction  errors  resulting  from  the  use  of  the 
regression  equations.   Of  course,  the  first-mentioned  type  causes  a 
decrease  of  resolution  in  the  prediction  equations.   The  prediction 
errors  have  already  been  discussed  to  some  extent  in  section  7. 

With  regard  to  data  acquisition,  all  data  recorded  for  this  study 
were  systematically  checked  for  errors  and  corrected  at  each  step  where 
transformations  or  transcriptions  occurred.   In  addition,  a  final  overall 
check  was  conducted  to  insure  that  all  accessible  errors  had  been  removed 
or  corrected. 

Centers  were  located  accurately  by  drawing  isobars  to  smaller 
intervals,  even  as  small  as  0.5-mb  intervals  where  necessary*   Neverthe- 
less, the  charts  were  small  in  their  physical  size,  and  errors  of 
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location  may  be  as  large  as  +  60  run. 

Some  error  exists  in  the  measurement  of  distances  in  this  problem 
since  all  distances,  including  the  distances  between  grid  points,  were 
measured  using  a  scale  assumed  true  at  latitude  45N.   This  is  the  central 
parallel  for  the  forecast  area. 

Since  the  charts  used  were  stereographic  projections  true  at  60N, 

the  map  factor,  m,  is 

m   =      1.86603  (6) 

1  -f-  sin  L 

If  one  assumes  a  movement  M  over  a  5-degree  latitude  segment  from  30  to 

60N,  use  of  the  mean-value  theorem  leads  to  the  percent  error 

-T7-    =     -   ~  (7) 

as  a  function  of  map  factor.   Errors  of  this  type  are  given  in  table  5. 
In  general,  the  scale  factor  contributed  only  a  small  percentage  error 
to  any  particular  computation  of  movement. 


Latitude               Map  factor  °L   Error 

30-35                  1.213  -11.0 

35-40                  1.160  -  6.1 

40-45                  1.114  -  1.9 

45-50                  1.074  1.7 

50-55                  1.040  4.8 

55-60                  1.011  7.5 

Table  5.   Errors  incurred  in  measuring  distance,  based  on 
5-degree  latitude  movements. 

Although  the  "t"  test  mentioned  near  the  end  of  section  7  indicates 
that  the  means  could  be  from  two  samples  drawn  from  the  same  population, 

this  does  not  preclude  the  possibility  that  the  independent  sample  might 

conceivably  be  from  a  different  population.  Possibly  the  most  plausible 
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explanation  is  given  by  Panofsky  and  Brier  [1958]: 

The  equation  with  the  most  predictors  will  not  necessarily 
yield  the  best  fit  to  the  second  sample.   The  reason  is  that 
the  longest  equation  may  have  actually  overfitted  the  first 
sample  and  ascribed  some  of  the  variation  due  to  small-scale 
fluctuations  to  one  of  the  predictors  "by  accident." 

In  light  of  the  above,  each  equation  of  the  intermediate  steps  of 
regression  is  shown  in  table  7.   In  this  manner,  future  research  in  the 
field  may  show  better  results  with  fewer  variables  as  suggested  by 
Panofsky.   One  can  calculate  the  significance  of  the  truncated  regression 
equations  by  using  the  appropriate  P.R.  variance  as  listed  in  table  3 
with  the  appropriate  degrees  of  freedom. 

Certainly  some  inaccuracy  results  from  the  fact  that  no  stratifica- 
tion of  the  variables  was  considered.   Thus,  for  example,  it  was  noted  in 
this  study  that  cold  and  warm  anticyclones  have  differing  movement 
characteristics.   It  would  appear  that  at  least  a  bimodal  distribution 
with  respect  to  mean- temperature  distribution  would  have  been  possible. 
9.   Conclusions. 

The  prediction  equations  are  at  least  an  improvement  over  known  exist- 
ing techniques.   They  seem  to  have  wide  applicability  over  a  vast  geo- 
graphic area.   The  past  history  of  an  anticyclone  track  is  not  required, 
nor  are  the  equations  limited  to  a  particular  flow  pattern  or  thermal 
distribution  aloft. 

It  is  felt,  however,  that  a  stratification  of  the  data  based  on  the 
thermal  character  would  have  resulted  in  a  greater  predictive  value  in 
the  final  equations.   For  example,  had  the  dependent  sample  been 
restricted  to  cold  surface  highs,  greater  predictability  would,  no  doubt, 
have  resulted.   It  is  felt  that  decreasing  the  southerly  extent  of  the 
forecast  area  to  40N,  thereby  partially  avoiding  the  subtropical  high 
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pressure  belt,  v>/ould  also  have  increased  the  predictive  value  of  the 
forecast  equations. 

The  results  suggest  that  further  investigation  along  these  lines 
should  prove  quite  rewarding. 
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APPENDIX  A 


TABULATION  OF  INDEPENDENT  PREDICTIONS 


Table  6  contains  predicted  values  of  N,  E  and  P  for  each  of  the  50 
cases  in  the  independent  data  sample.   In  addition,  verifying  values  are 
listed  as  well  as  the  forecast  errors  (N?  -  N) ,  etc.,  which  are  derived 
directly  from  the  computer  output.   Similar  remarks  apply  to  E1  and  E, 
and  the  component  error  is  (E'  -  E) .   The  vector  error  in  the  forecast 
is  equal  to  the  vector  sum  of  the  northward  and  eastward  component  errors, 
A  diagram  (fig.  6)  depicts  each  one  of  these  vector  errors. 

Negative  values  of  P  or  P'  represent  observed  or  forecast  anti- 
cyclone weakening.   The  error  is  equal  to  (P'  -  P) . 

All  displacements  are  measured  relative  to  the  center  of  the  basic 
square  grid  rather  than  the  projection  of  the  charts  from  which  the  data 
are  obtained. 
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APPENDIX  B 

SUMMARY  OF  FORECAST  PROCEDURE 

The  percentage  of  the  total  variance  of  a  predictand  explained  by 
each  selected  predictor  is  given  in  table  3.   In  general,  the  first  two 
or  three  predictors  explain  a  considerable  amount  of  the  variability; 
subsequent  predictors  explain  lesser  amounts. 

In  light  of  the  above,  where  time  is  the  primary  consideration,  the 
user  may  be  willing  to  sacrifice  some  of  the  predictive  value  of  the 
derived  equations  in  order  to  reduce  the  number  of  calculations  required, 
The  coefficients  for  the  truncated  equations  obtained  from  the  inter- 
mediate steps  of  regression  are  summarized  below  in  table  7. 


Table  7.   Summary  of  the  truncated  as  well  as  of  the  final  form  of  the 
regression  equations  for  northward  and  eastward  movement  and 
pressure  change  of  sea-level  anticyclones. 


N' 

E' 

P' 

k  =  3 

Const 

1225.88867 

Const 

-80.84535 

Const 

247.51972 

Z34 

-74 
T 
56 

-3.40210 

2.06476 

13.80024 

Z52 

?65 
^Z25 

-2.35890 

2.72546 

-1.98448 

P65 

A52 
Az45 

-.23060 

-.,01609 

.06027 

k  =  4 

Const 

-7.07663 

Const 

-844.11490 

Const 

356.69707 

Z34 
^74 

-2.34752 
2.30506 

A 
Z52 

6.72400 
-2.53248 

P25 
P65 

-.12058 
-.21709 

T 

-43 

-9.59201 

Z65 

3.07734 

Z52 

-.01637 

T 

X56 

15,14506 

Az25 

-2.83424 

AZ45 

.05371 

k  =  5 

Const 

6030.72598 

Const 

-1058.52912 

Const 

376.14332 

P22 
Z34 

|43 
56 

-6.30295 

* 

8.12394 

P25 
Az45 

-.11976 

-2.13926 

P63 
Z52 

Az25 

6.97667 

-.23663 

2.45926 

-2.32633 

.08134 

-13.74407 
17.32802 

2.91559 
-2.75370 

-.01671 
.04914 
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Table  7  (continued) 


N' 

E* 

P' 

k    =  6 

Const 

14646.77566 

Const 

5242.92812 

Const 

358.52542 

P22 

P97 
Z34 
f74 

?43 
56 

-5.42540 

A 

8.89200 

P25 

P65 

AP72 

A52 

^Z45 

-.12976 

-8.14462 

P84 

^•P63 

Z52 

A65 
Az25 

-6.50886 

-.21104 

-2.09395 

7.64845 

.14006 

2.14122 

-15.83553 

22.12236 

-2.22808 

3.10035 

-2.35194 

.08713 

-.01443 

.04065 

k   =  7 

Const 

11578.46055 

Const 

1537.06168 

Const 

333.43721 

P22 

P85 
P97 
Z34 
§74 
|43 
56 

-6.85958 

A 

8.83391 

P25 
P65 

Ap43 
Ap72 

t¥ 

7Z45 

Az92 

-.12157 

5.55932 

P84 

A85 

^p63 

Z52 

f65 
&z25 

-12.23171 

-.19594 

-9.84316 

-2.27834 

1.84021 

-14.81254 

23.40104 

9.41494 
7.2116 

-2.07604 
2.89902 

-2.44941 

.14064 
.09826 
-.01254 
.04010 
.01895 

k   =  8 

Const 

5060.13022 

Const 

320.55524 

A 

9.65318 

P25 

P65 

£P43 

ZiP72 

Z52 

f45 
Az92 

-.12599 

P33 
P84 

Ap63 

Z52 

az25 

-4.96103 

-.18277 

-12.37582 
10.98861 

.12548 
.10397 

7.03236 

-.02079 

-2.10596 

.01295 

2.93012 

.03841 

-2.46111 

.02271 

k   =  9 

Const 

4655.21111 

Const 

343.12419 

A 

10.44113 

P25 

^P43 

^P72 

AP76 

Z52 

Z74 

■    *  7 

A45 
^Z92 

-.12004 

P33 
P74 
p84 

A85 
Ap63 

Z52 
z ,_ 
A65 

Az25 

-7.80116 

-.21145 

7.60599 

.13452 

-16.25978 
10.44469 

.10122 
.07752 

6.32668 

-.02031 

-2.26007 

.01336 

3.03432 

.04002 

-2.60004 

.02344 

It  should  be  remembered  that  when  using  the  forecast  procedure, 
northerly  and  easterly  displacements  are  measured  relative  to  the  central 
point  of  the  square  grid  shown  in  figure  2.   Thus,  on  a  polar 


stereographic  projection,  easterly  displacements  are  not  measured  along 
latitude  circles,  but  normal  to  a  meridian  passing  through  the  anti- 
cyclone center. 

The  accompanying  worksheets  were  prepared  to  facilitate  the  employ- 
ment of  the  derived  prediction  equations.  In  addition,  they  are  used  to 
demonstrate  the  actual  procedure  followed  in  making  an  operational  fore- 
cast with  the  forecast  scheme.  The  case  in  question  is  one  included  in 
the  independent  sample.  At  1200Z  HJan59  this  sea-level  anticyclone  was 
located  at  40 . 3N  and  85. 7W. 

The  symbols  used  in  the  regression  equations  are  defined  below. 
Specific  values  of  predictors  refer  to  the  example  forecast.   Observation 
time  is  1200Z  HJan59.   All  subscripts  refer  to  grid-point  locations  as 
specified  in  figure  2. 

p_9       Surface  pressure  in  whole  millibars. 

/^  p, _       The  past  24-hour  surface  pressure  change  in  whole  milli- 
bars. 

z  ,       500-mb  height  in  tens  of  feet  with  the  ten-thousands  digit 
omitted  (e.g.,  18,770  is  recorded  as  877). 

Az/<-       The  past  24-hour  500-mb  height  change  in  tens  of  feet 
(e.g.,  18,670  to  18,800  is  recorded  as  +13). 

T,        Mean  temperature  in  whole  degrees  centigrade  in  layer  from 
surface  to  500  mb. 

A        The  longitude  of  the  surface  anticyclone  center  recorded 
to  the  nearest  tenth  of  a  degree. 

From  these  worksheets  the  24-hour  predicted  displacement  is  calcu- 
lated to  be  163  nautical  miles  east  and  462  nautical  miles  south.   This 
displacement  is  then  measured  along  the  square  grid  from  the  surface 
anticyclone  center  at  1200Z  HJan59.   The  actual  movement  was  190  nautical 
miles  east  and  530  nautical  miles  south.   The  forecast  change  in  central 
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pressure  was  -1  millibar.   Actually  no  change  occurred;  the  central 
pressure  remained  at  1033  mb. 

The  verifying  sea-level  pressure  chart  for  this  example  is  for  1200Z 
12Jan59. 
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NORTH  AMERICAN  ANTICYCLONE  PREDICTION  WORKSHEET 


SOUTH-NORTH  DISPLACEMENT 
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A  XI 
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Pi 
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constant   = 

FORECAST  S-N  DISPLACEMENT   : 
Direction  (  +  North,  -  South)   : 
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NORTH  AMERICAN  ANTICYCLONE  PREDICTION  WORKSHEET 
WEST -EAST  DISPLACEMENT 

X  ?£Z  x  10.44113     ,  ,'"M  iC'lM 

ft    ■' 
P33  A  -  ■ x        (-)7. 80116      = 

py4     __: x  7.60599  = 

Pg4                JC  /C x  (_)16. 25978  = 
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^P6       ^ x  6.32668  = 


z52  9  f'r x        (-)2. 26007      = 


zr.  ■'■''■'' x  3.03432      = 

65      
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Az25  "    '  x        (-)2, 60004  = 

constant  = 

FORECAST  W-E  DISPLACEMENT  : 

Direction  (  +  East,  -  West  )  : 
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NORTH  AMERICAN  ANTICYCLONE  PREDICTION  WORKSHEET 


CHANGE  IN  CENTRAL  PRESSURE 
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PROGRAM 

THE 

INSTITUTE 

OF  NAVIGATION 


TWENTY-FOURTH  ANNUAL 
MEETING 


19-21  JUNE,  1968 


MONTEREY,  CALIFORNIA 


Hosted  by 
THE  U.S.  NAVAL  POSTGRADUATE  SCHOOL 


THE  INSTITUTE  OF  NAVIGATION 

TWENTY-FOURTH  ANNUAL  MEETING 

MONTEREY  CALIFORNIA 

TUESDAY  -  18  JUNE,  1968 

REGISTRATION  -  Mark  Thomas  Inn  6:00  -  9:00  pm 

WELCOMING  RECEPTION 

Mark  Thomas  Inn  6:00  -  10:00  pm 

WEDNESDAY  -  19  JUNE,  1968 

SPEAKERS'  BREAKFAST  FOR  SPEAKERS  AND  SESSION 
CHAIRMEN  OF  THE  DAY  -  Mark  Thomas  Inn       7:15  am 

REGISTRATION  DESK  OPENS  KING  HALL, 

U.S.  NAVAL  POSTGRADUATE  SCHOOL  8:00  am 

TWENTY-FOURTH  ANNUAL  MEETING 

COMMENCES  9:00  am 

OPENING  REMARKS  BY  GENERAL  CHAIRMAN, 
W.  J.  TULL,   EXEC.  V.P.,   ION 

WELCOME  BY  SUPERINTENDENT,   U.S.  NAVAL 
POSTGRADUATE  SCHOOL,   REAR  ADMIRAL 
ROBERT  W.   McNITT,  USN. 

REMARKS  BY  ION  PRESIDENT,  DR  .  GENE  R.  MARNER 


I     —  SPACE  SESSION 

CHAIRMAN,   DR.  RICHARD  M.  HEAD, 
NASA,   ERC 


9:30  am 


Orbiting  Velocity  Meter  —  S.  R.   Sporn,   American 
Bosch  Arma  Corp.  ,    Garden  City,   N.Y. 

Irradiation  &  Manual  Navigation  —  Richard  F. 
Haines  and  William  H.  Allen,  NASA,   Ames, 
Moffett  Held,   Calif. 


COFFEE  BREAK 


10:30  -  10:45  am 


Sight  Reduction  Tables  for  Orbital  Navigation 
William  Devereux,    Kollsman  Inst.  Corp.  , 
Syosset,  N.Y. 

Manual  Onboard  Computation  Procedures  & 
Devices  for  Orbital  Navigation  &  Guidance  — 
Samuel  P.   Altman,   General  Electric  Co.  , 
Philadelphia,   Pa. 


WEDNESDAY  -  19  JUNE  (continued) 

LUNCHEON  —   Mark  Thomas  Inn  12:00  -  1:30  pm 

Speaker:  Jon  Lindbergh,   Manager 

Pac.  N.W.   and  Japan, 
Oceans  Systems,   Inc. 


II  -  MARINE  SESSION 

CHAIRMAN,  DEAN  CARLE.  MENNEKEN, 
U.S.  NAVAL  P. G.  SCHOOL 


2:00  pm 


A  Bench  Mark  Navigation  System  for  Surface  and 
Subsurface  Vessels  —  Gilbert  Fain,   Raytheon  Co. , 
Portsmouth,   R.I. 

Yacht  Navigation  —  Ben  Warriner,  Fort  Wayne,  Ind. 

Channel  Navigation  &  Docking  of  Supertankers  — 
Jack  Kritz  and  Morton  J.  Howard,  Kollsman  Inst. 
Corp.,   Syosset,  N.Y. 


COFFEE  BREAK 


3:30  -  3:45  pm 


Marine  Collision  Avoidance:    Human  Factor  Con- 
siderations for  the  Development  and  Operation  of 
an  Effective  Merchant  Marine  Radar  —  Thomas 
D.  Mara,   General  Dynamics,   Groton,  Conn. 

Navigation  Systems  &  Insurance  —  Walter  B.  Potts, 
Marsh  &  McLennan,   Inc.,  New  York,   N.Y. 


THURSDAY  -  20  JUNE,  1968 

SPEAKERS'  BREAKFAST  FOR  SPEAKERS  AND 
SESSION  CHAIRMAN  OF  THE  DAY  — 

Mark  Thomas  Inn  7:00  am 


III  -  GENERAL  SESSION  8:30  am 

CHAIRMAN,  WINFIELD  H.  ARATA, 
NORTHROP  NORTRONICS 


The  National  Plan  for  Navigation  — 
Alton  B    Moody,  FAA,  Washington,  D.C. 

Application  of  Satellite  Navigation  Techniques 
to  Marine  &  Air  Navigation  —  Joseph  Chernof, 
ITT  Aerospace,   San  Fernando,  Calif. 

Expanded  Application  of  Navy  Navigation  Satellite 
System  in  Marine  &  Air  Navigation  —  E.  F. 
Gallagher,   Navy  Astronautics  Group,   Point  Mugu, 
Calif. 


COFFEE  BREAK 


10:00  -  10:15  am 


USAF  Navigation  Satellite  Program  — 
Maj.  F.  M.  Charette,   SAMSO,   USAF, 
Los  Angeles,  Calif. 

An  Optimal  Inertial/Doppler-Satellite  Navigation 
System  —  R.  G.  Brown  and  L.  L.  Hagerman, 
Iowa  State  University,   Ames,   Iowa 


LUNCHEON  -  Mark  Thomas  Inn 


12:00  -  1:30  pm 


Speaker:  Brig.  Gen.   Robert  A.  Duffy,  USAF 

Deputy  for  Re-entry  Systems, 
SAMSO 


IV.  -  GENERAL  SESSION 


2:00  pm 


Flight  Test  Performance  of  an  Airborne  Omega 
Equipment  Set  —  G.  Rex  Morin,  Lear  Siegler, 
Inc.,    Grand  Rapids,   Mich. 

Computerized  Airborne  Omega  Navigation  — 
D.  B.  Daniel,  Northrop  Nortronics,   Hawthorne, 
Calif. 

Omega  —  Inertial  Hybrid  Receiver  —  Winslow 
Palmer,  Consultant,  Babylon,   N.Y. 


COFFEE  BREAK 


3:30  -  3:45  pm 


Precise  Positioning  of  a  Ship  at  Sea  Utilizing 
VLF  Transmissions  —  Method  G  Results  — 
R.  Bernstein,   G.   N.  Ruppert,   and  C .  O.  Bowin, 
IBM/WHOI,   Gaithersburg,   Md. 

Time-Difference  Position-Determination  System 
for  Aerospace  &  Terrestrial  Applications  —  J.  E. 
Gaffney,  Jr.  ,  IBM,  Federal  Systems  Division, 
Gaithersburg,   Md. 


ANNUAL      BANQUET 

OFFICERS1  OPEN  MESS 
U.S.  NAVAL  P.G.  SCHOOL 


RECEPTION:    6:00  pm 


DINNER:    7:30  pm 


PRESENTATION  OF  ANNUAL  AWARDS  BY  DR. 
GENER.  MARNER,  PRESIDENT,  THE  INSTITUTE 
OF  NAVIGATION. 


FRIDAY  -  21  JUNE,  1968 

SPEAKERS'  BREAKFAST  FOR  SPEAKERS 
AND  SESSION  CHAIRMAN  OF  THE  DAY 
Mark  Thomas  Inn 

V  —  AIR  SESSION 


7:00  am 


8:30  am 


CHAIRMAN,   JOSEPH  A.  PARINI, 
LEAR  SIEGLER,   INC. 

Separation  Assurance  Over  the  North  Atlantic  — 
Nathaniel  Braverman  —  FAA,   Atlantic  City,  N.J. 

Commercial  Navigation  for  the  1970's  —  R.  H. 
Waldman,   International  Airline  Navigators 
Council,  Quebec,  Canada 

Trends  in  Computer  Technology  —  Glen  M. 
Harold,   Sperry  Rand  Corp.  ,   St.  Paul,   Minn. 


COFFEE  BREAK 


10:00  -  10:15  am 


Integrated  Inertial  Velocity-Aided  &  Position- 
Aided  Aircraft  Navigation  —  Dr.  Boris  Danik, 
General  Precision,   Inc.,  Wayne,   N.J. 

Area  Navigation  —  A  Major  Capability  Improve- 
ment for  Air  Traffic  Control  —  Glen  A.  Gilbert, 
Glen  A.  Gilbert  Associates,   Miami,  Fla. 


LUNCHEON  -  Mark  Thomas  Inn 


12:00  -  1:30  pm 


Speaker:  Willis  M.  Hawkins,  V.P.  Science  & 

Eng.  ,  Lockheed  Aircraft  Corp. 


VI  —  AIR  SESSION 


2:00  pm 


The  Integration  of  Air  Navigation  Systems  — 
John  Larsen,   Consultant,   Annapolis,    Md. 

A  Navigation  Satellite  System  for  Use  with 
Army  Tactical  Aircraft  —  Keith  D.  McDonald 
and  Laurin  G.  Fischer,  Communications  & 
Systems,   Inc.,   Falls  Church,   Va. 

Optimum  Integration  of  Aircraft  Navigation 
Systems  —  William  Zimmerman,  Dynamics 
Research  Corp. ,   Stoneham,  Mass. 


COFFEE  BREAK 


3:30  -  3:45  pm 


Electronic  Terminal  Guidance  for  All-Weather 
VTOL  Operations  —  David  W.  Young,   Lockheed 
Electronics  Co. ,   Plainfield,   N.J. 

CAINS:    Inertial  Navigation  for  Carrier  Aircraft 
in  the    1970's  —   R.    S.    Vaughn,    Naval  Air 
Development  Center,  Johnsville,  Pa. 

MEETING  ENDS  ABOUT  4:45  pm 


TWENTY-FOURTH  ANNUAL  MEETING 

General  Chairman W  .  J .  TULL 

Technical  Chairman LOREN  E.  DeGROOT 

Arrangements  Committee: 

Chairman RALPH  B .  NEAL 

Assistant  Chairman-Registration  . .  . 

...  LT.  COL.  PHILLIP  E.  HUKILL,   USA 

Assistant  Chairman WILLIAM  J.  PYESTNER 

Assistant  Chairman THOMAS  E.  DAWSON 

Assistant  Chairman-Publicity HAL  WALTON 

Ladies'  Program JOYCE  NEAL 

***** 


REGISTRATION 

A  registration  form  is  attached  (reverse  of  this  panel). 
Please  register  and  sign  up  for  meal  functions  early.  If 
the  Treasurer  is  so  requested  not  later  than  48  hours  in 
advance  of  the  meeting,  payment  for  registration  fee 
and  meals  will  be  refunded. 


HOTEL  ACCOMMODATIONS 

ION  Headquarters  will  be  located  at  the  Mark  Thomas 
Inn,  P.O.Box  1711,  Monterey,  California  where  a  block 
of  rooms  is  being  held  for  The  ION  until  5  June.  MON- 
TEREY IS  A  POPULAR  RESORT  AND  GOLFING  AREA, 
SO  MAKE  YOUR  RESERVATIONS  EARLY.  Room  rates 
are  as  follows: 


The  Institute  of  Navigation 

NATIONAL  OFFICE:     711   14th  Street,   N.W.,   Suite  912 
Washington,  D.C.   20U05 

PURPOSE  OF  THE  ION 

The  Institute  of  Navigation  is  a  scientific,  non-profit 
organization,  founded  in  1945.  Its  programs  are  directed 
toward  elevating  standards  of  navigation  by  coordinating 
the  knowledge  and  achievements  of  practicing  naviga- 
tors, scientists  and  those  involved  in  the  development 
and  production  of  navigational  equipment. 


Single 
Twin 


-  $  9.00  to  $18.00 

-  $13.00  to  $22.00 


When  requesting  a  reservation,  state  that  you  will  be 
attending  The  ION  Annual  Meeting. 


FUTURE  ION  MEETINGS  SCHEDULED 

6-7  November,    1968 

NATIONAL  MARINE  NAVIGATION  MEETING 
New  Orleans,  Louisiana 


President DR .  GENE  R .  MARNER 

Executive  Vice  President WILLIAM  J . .  TULL 

Treasurer COL.  WILL  O.  BRIMBERRY,  USAF 

General  Counsel ROBERT  REED  GRAY 

Eastern  Regional  Vice  President E.  S .  KEATS 

Central  Regional  Vice  President... 

...  MAJ.  WM.  L.  POLHEMUS,  USAF  (Ret.) 
Western  Regional  Vice  President WINFIELD  ARATA 


Editor 


THE  JOURNAL 
FREDERICK  FRANKLIN 


15-16  January,    1969 

FIFTH  UNCONVENTIONAL  INERTIAL  SENSORS 
SYMPOSIUM  (Classified)   —   NAVAL   APPLIED 
SCIENCE  LABORATORY,    BROOKLYN,    N.    Y. 
CO-SPONSORED  WITH  THE  NAVY   AND   AIR 
FORCE. 

5-6  March,    1969 

NATIONAL  SPACE  MEETING 
Houston,  Texas  (Tentative) 


Executive  Director  . .  . 

...  CAPTAIN  ROSS  E.  FREEMAN,  USN(Ret.) 


r 


